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Gaseous drag on a gravitational perturber in Modified 
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ABSTRACT 

We calculate the structure of a wake generated by, and the dynamical friction force 
on, a gravitational perturber travelling through a gaseous medium of uniform density 
and constant background acceleration g ext , in the context of Modified Newtonian 
Dynamics (MOND). The wake is described as a linear superposition of two terms. The 
dominant part displays the same structure as the wake generated in Newtonian gravity 
scaled up by a factor /i _1 (<7ext/»o), where ao is the constant MOND acceleration and 
/i the interpolating function. The structure of the second term depends greatly on the 
angle between g cxt and and the velocity of the perturber. We evaluate the dynamical 
drag force numerically and compare our MOND results with the Newtonian case. 
We mention the relevance of our calculations to orbit evolution of globular clusters 
and satellites in a gaseous proto-galaxy. Potential differences in the X-ray emission of 
gravitational galactic wakes in MOND and in Newtonian gravity with a dark halo are 
highlighted. 

Key words: hydrodynamics - galaxies: haloes - galaxies: interactions - galaxies: 
kinematics and dynamics - galaxies: structure - gravitation 
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1 INTRODUCTION 

According to the Brownian theory, a 'macroscopic' parti- 
cle moving in a fluid experiences a fluctuating force as a 
consequence of the graininess of the fluid, resulting in dy- 
namical friction (DF) and diffusion in the velocity-space. 
Understanding the nature of the DF force experienced by a 
gravitational object that moves against a mass density back- 
ground is of great importance for describing the evolution of 
gravitational systems and the exchange of angular and lin- 
ear momentum between gravitational subsystems. There are 
two limits of interest in astrophysics: pure collisionless media 
and fully collisional gas. If the pertuber only interacts with 
the surrounding particles gravitationally, either the medium 
is collisionless or collisional, the DF force is the result of the 
formation of an overdensity wake behind the perturber. 

The DF depends on the nature of the force between the 
perturber and the background particles. Departures from 
Newtonian gravity have been proposed as an alternative to 
dark matter in galaxies (e.g., Modified Newtonian Dynam- 
ics; Milgrom 1983). Ciotti & Binney (2004) found that in 
collisionless backgrounds, the characteristic timescale for DF 
in Modified Newtonian Dynamics (MOND) is much shorter 
than in the Newtonian equivalent system with dark mat- 
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ter. If confirmed, this result may have strong astrophysical 
implications and may be used to distinguish between dark 
matter and modified gravity. Using the MOND scaling found 
in Ciotti & Binney (2004), the globular clusters in Fornax 
dwarf spheroidal galaxy should have spiraled into the centre 
of Fornax in ~ 0.5 Gyr, forming a visible galactic nucleus 
(Sanchez-Salcedo et al. 2006). 

Ciotti & Binney (2004) derived the MOND DF 
timescale in a restrictive plane-parallel geometry, which is 
a too specialized case, and in a non-standard statistical for- 
mulation. These caveats led Nipoti et al. (2008) to carry out 
fully N-body simulations of the dynamics of bars in MOND 
systems. Their simulations confirmed the scaling of the DF 
timescale predicted in Ciotti & Binney (2004) . Still, it is not 
entirely clear how the background responds to the pertur- 
bation, and the description of the structure of the density 
wake induced by a small and massive perturber is lacking. 
Since MOND is a nonlinear theory, the gravitational accel- 
eration induced by the perturber is affected by the presence 
of an external gravitational field. In particular, the DF force 
is expected to depend on the angle between the direction of 
the mean field and the velocity of the perturber. Our aim is 
to understand the nature of DF in MOND. 

In this paper we consider gaseous media and study the 
structure of the density wake induced by a gravitational per- 
turber in MOND. Gaseous DF has found many applications 
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in astrophysics, ranging from protoplanet accretion all the 
way to the motion of galaxies in clusters (e.g., Ostriker 1999; 
Kim 2007; Conroy & Ostriker 2008) . In a pure baryonic uni- 
verse as in MOND theory, the relative role of gas becomes 
even more important because all the contribution to the 
mass is accounted for by the gas and the stars. Moreover, 
the complementary view of hydrodynamics may provide new 
insights into the differences and analogies between Newto- 
nian and modified gravities. Although the details on the 
internal density structure of the induced wake may depend 
on whether the system is fully collisional or collisionless, the 
fluid limit is useful for understanding the main conceptual 
features introduced by the change of the gravity law. 

The linear response of the gaseous medium to a gravi- 
tational pcrturber in Newtonian gravity is well-documented 
(Dokuchaev 1964; Ruderman & Spiegel 1971; Just & Kegel 
1990; Ostriker 1999; Sanchez- Salcedo & Brandenburg 1999, 
2001; Kim & Kim 2007; Kim et al. 2008). In all these works, 
a minimum radius r m i n in the Coulomb logarithm was intro- 
duced in order to regularise the gravitational potential of a 
point mass. Although the formulae were derived for rectilin- 
ear orbits in homogeneous and infinity media, simple 'local' 
extensions have been proven very successful in more realistic 
situations, e.g., smoothly decaying density backgrounds or 
when the perturber is moving on a circular orbit (Sanchez- 
Salcedo & Brandenburg 2001; Kim & Kim 2007). 

The paper is organized as follows. In §2, we discuss the 
basic concepts on the ideal problem of a point particle mov- 
ing at constant speed within a uniform gas in MOND (the 
Bondi-Hoyle problem). In §3, we outline the linear deriva- 
tion of the basic equations for calculating the steady-state 
density wake generated by an extended body. In §4 we de- 
scribe the structure of the resulting wake and evaluate the 
DF force exerted on the perturber. We then discuss some 
implications of our results in §5. Finally, we conclude in §6. 



the MOND acceleration g is related to the Newtonian ac- 
celeration, g N by the algebraic relation p(\g\/a )g = g N 
(Brada & Milgrom 1995). At variance with the Poisson equa- 
tion, MOND is a non-linear theory. This implies that the 
gravitational potential generated by the perturber depends 
in a complex way on the external field g ext in which it is 
immersed. The external field is to be thought of as the field 
of an enveloping system in the absence of the perturber. 
For example, the perturber can be a galaxy in the external 
field of a cluster of galaxies, or a dwarf spheroidal galaxy in 
the field of its parent galaxy. In MOND, it is crucial to in- 
clude the external field g cxt to describe DF on a perturber. 
The case g oxt = is not very relevant when considering DF 
because it corresponds to a situation where all the surround- 
ing gas is bound to the perturber, whereas DF concerns the 
interaction with unbound distant particles having large im- 
pact parameters. MOND permits to have unbound particles 
if g ext 7^ 0. For these reasons, we will assume that g oxt is a 
no-null vector field. 

Along this paper we will assume that the unperturbed 
gas density is constant all over the space and that the body 
moves in a uniform rectilinear orbit, which seems to be in 
contradiction with the inherent assumption that there is an 
external gravitational field. Note, however, that this approx- 
imation is amply used even in Newtonian dynamics where in 
most of the cases (e.g., a galaxy inside the galaxy cluster) the 
gravitational body is immersed in a stratified medium. The 
'local' approximation, that is, estimating the drag force at 
the present location of the perturber as if the medium were 
homogeneous but taking appropriately the Coulomb loga- 
rithm, has been proven very successful in both collisionless 
and collisional systems (e.g., Sanchez- Salcedo & Branden- 
burg 2001, and references therein). So far we are not inter- 
ested in including the complications of density gradients in 
the unperturbed surrounding medium. 



2 THE BONDY-HOYLE PROBLEM IN MOND 

We consider a gravitational point particle at the origin of 
our coordinate system, surrounded by a gas whose velocity 
far from the particle is 



= Mc x z, 



(1) 



where c x is the sound speed of the gas at infinity and M is 
the Mach number. We are interested in the wake produced 
by the gravitational interaction with the perturber in the 
context of MOND. 

Bekenstein & Milgrom (1984) suggested a Lagrangian 
theory and also a nonlinear differential equation for the non- 
relativistic MOND gravitational potential produced by a 
mass density distribution p: 



p 



fflo 



4ttGp, 



(2) 



where a ~ 10~ 8 cm s~ 2 is an universal acceleration and 
p(x) is a monotonic and continuous function with the prop- 
erty that p(x) ~ x for i « 1 (deep MOND regime) and 
p{x) ~ 1 for x >• 1 (Newtonian regime). Only for very spe- 
cial configurations (one- dimensional symmetry -spherical, 
cylindrical or plane symmetric systems- or Kuzmin discs) 



2.1 The far-field equation: external dominated 
field 

Far away enough from the object, the change in the exter- 
nal gravitational potential can be treated as a linear per- 
turbation. This is valid as soon as the physical size of the 
perturber is supposed to be small as compared to the charac- 
teristic length-scale of the medium. Using the subscript for 
referring to unperturbed quantities, the linear field equation 



[a*o (l + £<)£<)) ■ V$i] = AnG Pl 



(3) 



where <3?i is the change in <£> produced by an increment pi 
in p (Milgrom 1986). Here po = Mffext/ao), Lo is the log- 
arithmic derivative of po (in the unperturbed system) and 
So represents a 3 x 3 matrix with elements £o,i,j = eo,ieoj, 
where eo is the unit vector in the direction of g cxt , which is 
r dependent. As Milgrom (1986) pointed out, $i satisfies an 
equation analogous to the electrostatic field equation with 
an inhomogeneous and anisotropic dielectric tensor. 

Since we are not interested in the distortions in the 
wake by tidal forces, we will assume that g oxt is a constant 
vector all over the space. This assumption also facilitates a 
comparison with the more familiar Newtonian case. If so, po, 
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Lo and eo are also constant and the linear field equation for 
a constant external field can be written as (Milgrom 1986): 

It can be seen that the potential becomes Newtonian 
(but with a larger effective gravitational constant) and 
anisotropic. We must note that, although the above equa- 
tion was derived for |V3>i| <§C <? C xt, it is also valid when 
ffoxt ^ cto, regardless the value of | V$i|. In this limit p,o ~ 1 
and Lq « 0; the conventional Poisson equation is recovered. 



2.2 The Bondi-Hoyle radius 

In the classical problem of a slow Brownian particle in a 
fluid, the field particles are assumed to form a heat bath. 
That is, they all stay close to thermal equilibrium, despite 
the presence of the Brownian particle and, hence, the equa- 
tion of motion is solved by perturbation. In the case of a 
point gravitational perturber immersed in a perfect gaseous 
medium, the Bondi-Hoyle radius defines the region where 
the response of the gas is linear. In Newtonian dynamics the 
Bondi-Hoyle radius is r B H = GM/(4(1 + M 2 )). Stream- 
lines whose impact parameter is less than 2tbh, will bend 
significantly and pass through a shock. Hence, within Itbh 
it is not any longer a small perturbation. In order to regu- 
larise the gravitational potential of a point mass, one has to 
introduce a minimum radius in the formulae of the DF drag 
(r min ~ 2r Bff ). 

In the Appendix A, we estimate the Bondi-Hoyle radius 
for a point mass in MOND. It is shown that the MOND 
Bondi-Hoyle radius is larger than in Newtonian gravity by a 
factor between p^ 1 and A t ( J 1 (l + Lo) 1 ^ 2 , depending on the 
angle between the velocity of the particle and the external 
field. To get a sense of values of tbh in typical cases, consider 
a galaxy of 5 x 10 11 M Q orbiting supersonically M ~ 1.5 in a 
cluster of galaxies with a sound speed of intracluster gas of ~ 
1500 km s _1 . The Bondi-Hoyle radius in this case is tbh 
/!q 1 kpc. In a typical galactic cluster po = 0.3-1 (see, e.g., 
fig. 7 in Sanders & McGaugh 2002), hence r B H < 1-3 kpc, 
implying that, if the interaction with the intracluster gas is 
merely gravitational, the linear approximation is satisfactory 
for studying the gaseous wake even quite close to the galaxy. 
For extended perturbers with characteristic size much larger 
than 2tbh, the flow is essentially laminar at any location. In 
the remainder of the paper we will describe the perturbation 
on the gas using linear theory. 



where the internal acceleration gi nt ~ GM p /(porp) >C ffcxt 1 - 
The following density-potential pair, which corresponds to a 
"modified" Plummer model, is an exact solution of Eq. (4) 



ft>0) = 
Mr) = 



3Mri 



AnVT+Lo- (x* + + _£i_ + r 2j 

1 GM_ 

£to ^/(l + L Q )(x 2 +y 2 ) + z 2 + ff 



5/2 ' 



(5) 



(6) 



where r p and r p 



(1 + LoY^r-p are the characteristic radii. 
Note that for Lo = and po = 1 (Newtonian limit) , it corre- 
sponds to the classical spherical Plummer model. Sometimes 
it is useful to express p p in terms of the central density of 
the object, p c , as follows 



P P (r) = P, 
where 

Pc = 



1 + ^7 + ^ + 



3M 



-i -5/2 



(l + io)r 2 



47r(l + L ) 1 /2 r 3- 



(7) 



(8) 



The central density can be written in terms of the central 
density for the spherical Plummer model in classical New- 
tonian gravity, p c , as p c = p c /(l + io) 1/2 . 

The selection of a Plummer model was for analytical 
purposes. As long as the size of the system is much larger 
than the perturber, the structure of the wake and the drag 
force experienced by the perturber are not expected to be 
sensitive to the details of the potential close to the body. 

In order to approach this problem analytically, we study 
the simplest case: an externally dominated perturber. In real 
life, there are some Galactic dwarf spheroidal galaxies that 
are known to be in this regime (e.g., Milgrom 1995; Sanchez- 
Salcedo & Hernandez 2007). It is likely that some subclusters 
and groups of galaxies, with low internal accelerations, em- 
bedded in a main massive galaxy cluster (such as the NGC 
4911 group in the Coma Cluster) lie also in this regime (e.g., 
Sanders & McGaugh 2002). 

In Eqs (5) and (6), the external field was taken along 
z and, therefore, in the same direction as the incident flow 
(see Eq. 1). In this section we will focus on this axisym- 
metric case. The derivation of the equations when g oxt is 
perpendicular to is postponed up to §4.1.3. These two 
situations brackets a general case where the external field 
has an arbitrary angle with respect to the velocity of the 
flow at infinity. 



3 FORMULATION 

3.1 Modelling the perturber 

Our aim is to study the large-scale gravitational perturba- 
tion induced by a small perturber travelling through a much 
larger system. We obtain that, beyond a certain distance 
from the perturber, |V3>i| <iC <? e xt- In other words, the far- 
field wake is expected to be in the external field-dominated 
regime. In order to highlight the differences between gen- 
uine MOND and Newtonian gravity, we restrict our consid- 
erations to situations in which the potential is dominated 
by the external field everywhere, i.e., we will consider an 
extended perturber of mass M and characteristic size r p , 



3.2 Linear equations in an external dominated 
field 

In the following, we give the linear derivation of the wake in 
a medium with unperturbed density po and adiabatic sound 
speed Coo, ignoring gas self-gravity and any magnetic fields. 
As stated in Eq. (1), it is assumed that the gravitational 
perturber is seated at the origin of our coordinate system 
and the gas velocity far from the perturber is Woo = McooZ. 
In the axisymmetric case, B ro and g ext are parallel. We are 

1 For a point-like perturber, one can always find a vicinity of the 
body where the inequality g lnt <C <7cxt is not achieved. 
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interested in the steady-state density enhancement T>(r) = 
(p — po)/po produced by the gravitational interaction with 
the perturber. The steady-state linearized basic dynamical 
equations for adiabatic perturbations p = p + p' and v = 
Woo + v' are: 



p V«' + Woo ■ Vp' = 0, 



and 

(«oo 



V)w' = -^Vp' 

Po 



(9) 



(10) 



Our strategy is to eliminate v 1 everywhere. By substituting 
equation (9) in the divergence of equation (10), we obtain 
that T> satisfies the differential equation 



LV = ^-V 2 $„ 



(11) 



where L is the linear differential operator 

The operator L arises frequently in fluid dynamics (e.g., Lan- 
dau & Lifshitz 1959). 

In Eqs (9)- (12) we have not specified the law of gravity. 
In MOND, the Laplacian of the potential in Eq. (11) can 
be expressed in terms of p v using the MOND field equation 
(4). The equation for V becomes 

LD = -^ + ^5%. (13) 

AtocJo cio oz 2 

From now on, it will be convenient to use the fol- 
lowing dimensionless variables: x = x/r p , y = y/r p and 
z = z/(l + Lo) 1 / 2 r J) . The hat symbol over a certain variable 
X will be used to denote that \ is written with x, y and z 
as the arguments 2 , e.g., D(x, y, z) = D(x, y, z). The second- 
order derivative in Eq. (13) can be performed as soon as 
the potential is known. Evaluating the second-order deriva- 
tive of the potential in Eq. (13) using the potential given 
in Eq. (6), and rearranging and grouping the terms, the so- 
lution V can be expressed as a linear superposition of two 
contributions V — V\ + T>2- Each one satisfies the following 
differential equations: 



LDi = 
and 
LV 2 = 
with 
9i(r) 

32(f) 



4?rGp c r 2 
To] 2~^gi(r), 



pocj 



AirGp c r 
-To T -^g 2 (r), 

Pocto 



[1 + r 2 ] s /2 ' 
2z 2 - x 2 



where f 2 



To 



[1 + 

= x 2 



^215/2 



+ y 2 + z 2 , 



3(1 + Lo) ' 



(14) 



(15) 



(16) 



(17) 



(18) 



2 At this stage it is obvious that this transformation does not 
conserve mass in the sense that if p is a density field, then 



and the operator L is 
dx 2 dy 2 

with 

2 _ M 2 + Lo 



dz 2 ' 



(19) 



(20) 



According to Eq. (18), To varies from in the Newtonian 
regime (Lo = 0) to 1/6 in the deep MOND regime (i.e. Lo = 
1). From Eq. (20) it can be seen that M. e g > M in the 
subsonic case, M. e g < M. in the supersonic case, and M e s = 
M = 1 at the transonic velocity. 

We recall that the equation for V in the Newtonian case 

is 



LV = — 



4-rrGpc 



^ffi(f), 



(21) 



which is naturally recovered from the above equations just 
by taking po = 1 and Lo = (so that M e s = M). The com- 
ponent T>\ obeys a differential equation similar to the New- 
tonian case and hence is identical to the wake induced by a 
perturber with mass density pi = (l-T„)(l+L )" 1/2 p c gi(f ) 
in conventional Newtonian gravity, once G has been replaced 
by a larger effective value G/po- The profile p\ corresponds 
to the classical (spherical) Plummer model in the trans- 
formed coordinates (x, y, z), multiplied by a form factor be- 
tween 0.59 (deep-MOND limit) and 1 (Newtonian limit). In 
analogy to the Newtonian case (Eq. 21), a fictitious New- 
tonian perturber with a pseudo-density mass density distri- 
bution p2 = Topc<?2 generates a component identical to T>2- 
We refer to p2 as pseudo-density because it may take nega- 
tive values. Interestingly, the total mass associated with this 
distribution is 



M 2 =VV 



/oo foo rc 
-00 J — 00 J — < 



p2 dxdydz = 0. 



(22) 



Our natural choice was to adopt the Cauchy principal value 
of the integral, which we denote by the symbol W (e.g., 
Mathews & Walker 1970). Note that p2 decays more slowly 
with radius than pi . Hence, the pseudo-density p2 becomes 
larger than pi at large radii. The importance of the contri- 
bution of T>2 to the wake and to the gravitational drag is 
difficult to forsee without a quantitative study. 



3.3 Formal solution 

The solution of Eqs (14) and (15) can be obtained using 
the retarded Green's function, which may be derived follow- 
ing different paths (e.g., Just & Kegel 1990; Ostriker 1999). 
In particular, the Green function can be found after Fourier 
transforming and imposing causality when choosing the con- 
tour of integration in the complex plane for M e g > 1 (e.g., 
Just & Kegel 1990; Furlanetto & Loeb 2002). In the steady- 
state, the perturbed density fields T>\ and T>2 are 



f pdxdydz ^ f pdxdydz. 



V 1 

T>2 

where 



1 - Tp Gp c r; 
(l + Lo) 1 / 2 ^ cio J 

T Gp c r 2 
(l + Lo) 1 / 2 ^ 



(23) 
(24) 
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Ii(r) = / d 3 



7,(25) 



for i = 1, 2, with /3 2 ff = X 



1 and 



' 2 if Mrff > 1, and 

2 - z' + f3 cB [R 2 + R' 2 - 2i?i?'cos(9'] 1/2 > 0; 
1 if M e s < 1; 
otherwise. 



4 RESULTS 

4.1 The density structure of the wake 

4-1.1 The component T>\ in the axisymmetric case 

Figure 1 shows the integral Ji, which is proportional to 
the perturbed density D\, in the z-R plane, of a MeS = 
0.8, 1.13, 1.5 body. In the deep MOND regime, they corre- 
spond to physical Mach numbers of 0.53, 1.25 and 1.9, re- 
spectively. So far we are only interested in the 'far-field' per- 
turbed density, hence we will not delve into details regarding 
the near-field (within a few core radius from the perturber). 
A subsonic perturber generates a density distribution with 
contours of constant density corresponding to similar ellipses 
with eccentricity A^eff ■ For supersonic motions, however, the 
region of perturbed density is confined within the rear Mach 
cone, dragged by its apex by the perturber. The surfaces of 
constant density within the wake correspond to hyperbolae 
in the z-R plane, with eccentricity e = Alcft. This is ex- 
pected because, as we show in Section 3.2, the equation for 
£>i has the same form as in the Newtonian case with Al c ff , 
once the density distribution of the perturber is rescaled by 
a factor (1 — Tq)(1 + Lo)~ 1//2 , and G is replaced by G/fio- 
Using this analogy, we can take advantage of the analyti- 
cal results of Ostriker (1999) to find the 'far-field' perturbed 
density: 



23i 



To ii^CGhhrl/ci 



(l + Lo) 1 ^ 



z 2 - Pl s R? 
Ho 1 S"GM/(Z 



(26) 



(27) 



where £" = 2 for supersonic perturbers and 1 in the subsonic 
regime. In Eq. (27) we used that Mi = p c J gidxdydz = 
M/rp. From the equation above, we see that the isoden- 
sity contours are z 2 + R 2 (l — M 2 ) =const, i.e. ellipses or 
hyperbolae with eccentricity e = A4. In the Appendix B we 
reconsider the perturbation as a time-dependent rather than 
a steady state problem. 
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Figure 1. Distribution of the dimensionless integral I\, in log- 
arithmic scale, for Af e fl = 0.8, 1.13 and 1.5 in the axisymmetric 
case. The box size for A'loff = 1-13 is larger in order to have the 
same dynamical range in density. 



4-1.2 The component T>2 in the axisymmetric case 

Figure 2 shows I2 for the same three effective Mach num- 
bers as previously considered (Meg = 0.8, 1.13 and 1.5). In 
the subsonic regime, I2, with a bipolar structure, is positive 
along the z axis and negative in the perpendicular plane. For 
supersonic perturbers, an overdense bump at the head of the 
perturber is generated. Interestingly, I2 displays a drop in 
density along the surface of the Mach cone. This negative 
jump in the surface of the Mach cone is very remarkable at 
MeS = 1.13 and dilutes at larger effective Mach numbers. 



Our calculations show that at M c s = 1-75 there is no jump 
in the Mach cone surface. For perturbers with M c s > 1.75 
the density jump in the Mach cone becomes positive. We 
see that, in general, I2 may take values comparable to ii. 
We note, however, that the real density D2 is related to I2 
through a factor T < 1/6 (see Eq. 24). 

Figure 3 contains the superposition (1 — To)h + T0/2, 
which is proportional to T>, in the deep MOND regime 
(i.e. Lq = 1 and To = 1/6). In the subsonic regime, the 
inclusion of the component T>2 plays an important role. The 
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resulting contours of isodensity can be fitted by ellipsoids 
denned by the equation R 2 + z 2 /q 2 =const, where q is the 
flattening parameter of the distribution in the physical R-z 
plane. According to our discussion in the preceding section, 
for M c fi = 0.8 and Lo = 1 (i.e. M = 0.53), we know that the 
isodensity contours of T>\ have q = 0.85. The superposition 
of components L>i plus T>2 generates a density distribution 
with q — 1.23 (in the deep-MOND regime). In order to have 
an axis ratio of q = 0.85, the body should be moving at 
M = 0.85. For M = 0.2, the flattening parameter is 0.98 
if only T>\ is considered, and becomes 1.38 when T>2 is also 
taken into account. As a general conclusion, the wake for a 
subsonic body is more flattened along the direction of mo- 
tion in MOND than in Newtonian dynamics. 

In the supersonic deep-MOND case, the overdensity 
head is still visible when both components are added. The 
inclusion of T>2, however, does not change significantly the 
structure of the wake within the Mach cone for A4 c s < 1.5. 
At larger Mach number the contribution of T>2 relative to 
T)\ becomes less and less important. Since the structure of 
T>\ is a scaled version of the wake generated in the New- 
tonian case, the difference between the structure of a wake 
generated by a small highly-supersonic perturber in MOND, 
would be likely too subtle to be distinguished from a per- 
turber of fictitious mass (1 — To)^ 1 in Newtonian gravity. 
Differences only appear in the vicinity of the perturber. 



4-1.3 External field orthogonal to the velocity flow 

Suppose now that the external field is along the x-axis. By 
denoting now x = x/r p , y = y/r p , z = z/r p , and by rea- 
soning entirely analogous to that leading to Eqs (14)-(20), 
one finds that when the external field is perpendicular to 
the velocity flow, M e s, gi and §2 are given by 



M c s = M, 



9i{r) = 
and 

§2 (r) = 



1 + 



1 + L 



,-2,-2 

+ y 



-, -5/2 



2x 2 



1+Lo 



-2 

y 



9i- 



(28) 
(29) 

(30) 



In this case, the dependence on Lo does not factorize so 
that one needs to calculate the integrals (25) for each pair 
(M cB ,L ). 

We found numerically that T>\ is very axisymmetric 
around the x axis; the effect of the gravitational dilation 
on T>i is small. For instance, when adopting Lo = 1, the an- 
gular variations of X>i are less than 2%, 2%, 8% and 13% for 
M = 0.4, M = 0.8, M = 1.25 and M = 1.8, respectively. 
T>\ is almost undistinguishable (differences of ~ 5 percent or 
less) from its counterpart in the axisymmetric case, and thus 
they are not shown. In particular, the opening angle of the 
Mach cone for supersonic perturbers in physical coordinates 
(z,R), is the same as in the axisymmetric case. 

Unlike T>\ , T>2 is expected to depart from axisymmetry 
about the z axis. We will focus again on the deep-MOND 
limit (Lo = 1). Figure 4 shows I2, which is proportional 
to T>2, in three perpendicular planes: y = 0, x = and 
z = 10r p , for M = 0.53 and M = 1.25. When the perturber 
moves subsonically, the structure of T>2 in the y = plane 
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Figure 2. Profile of the dimensionless integral I2, in the axisym- 
metric case. The scale is linear scale. The effective Mach numbers 
are the same as in Fig. 1. 



looks pretty much like in the axisymmetric case after a rota- 
tion of 7r/2 (compare Fig. 2 and Fig. 4). In the x = plane, 
however, the density map is notoriously different. It clearly 
shows that the configuration T>2 has not an axial symmetry 
around the z axis. 

For subsonic perturbers, T>2 may be able to change the 
flattening parameter of the wake. As an example, consider 
A4 — 0.8 and Lo = 1. In this situation, T)\ has isodensity 
contours with q = 0.6. If the contribution of T>2 is added, 
the composed wake exhibits q — 0.43 in the y — plane, 
and q — 0.54 in the x — plane (not shown). 
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Figure 3. Density perturbation profiles, p/po, in units of 
Pq 1 [GpcTp/V^c^o] in the axisymmetric case and with Lq = 1. 
The scale is logarithmic. 



For a supersonic perturber, the overdensity regions in 
T>2 are not located any longer at the head of the body; re- 
gions at the front exhibit a decrease in density. Zones with 
density depletion, that is T>2 < 0, can be found at both 
downstream and upstream. The maximum density enhance- 
ment in 2?2 appears smaller, by a factor ~ 4 at M = 1.25, 
than in the axisymmetric configuration. Figure 4 shows the 
complexity of the structure in the (x,y) plane (lower pan- 
els). The isodensity contours in that plane turn up to be 
elongated along the a;-axis. 



4.2 The gravitational drag on the body 

Once the structure of the wake T> is constructed, it is 
straightforward to evaluate the drag force exerted on the 
perturber by its wake. By symmetry, the only non- vanishing 
component lies along the z axis. In particular, the drag force 
in the axisymmetric case is 

GMp (' (' , .„„ Vz 



F DF = 2tt- 



IM> 



dzdRR 



((l + L )i? 2 + z 2 +r2)3/2 



(31) 



If Fdf,i and Fdf,2 denote the contribution to the drag force 
by components T>\ and T>2, respectively, we have: 



Fdi 



Fdfa + Fdf,: 



Fdfa = — 



To G 2 M 2 P0 



2p 2 (l + L ) c% 



Fdi 



where 



To 



G 2 M 2 p 



2p 2 (l + L ) 



T 2 



Ti = 



// 



dzdR- 



FRz 



(32) 
(33) 

(34) 
(35) 



(l + R 2 + i 2 ) 3 / 2 

We first estimate the relative contribution of Fdf,2 as 
compared to Fdf,i- For supersonic perturbers, Ti and T\ 
were calculated by carrying out the integration in Eq. (35) 
over all our box domain. Hence, the i-integral has up- 
per/lower limits z = ±25 and the R- integral has limits 
R — 0, 25. This implies that our steady-state wakes have 
an extent along z of z max ~ 25r p to 35r p , depending on Lq, 
and a Coulomb logarithm In A ~ ln« ma x/2r p ~ 3 (Sanchez- 
Salcedo & Brandenburg 1999). 

For a purely steady-state density perturbation, the 
forward-backward symmetry in the subsonic case argues 
that zero net force acts on the perturber. However, Ostriker 
(1999) noticed that this conclusion is misleading because, 
although complete ellipsoids exert no net force on the per- 
turber, there are always cut-off ones within the sonic sphere 
that exert a gravitational drag. Thus, we estimate the gravi- 
tational drag on a subsonic perturber by integrating Eq. (35) 
over the largest sonic sphere contained in our computational 
domain. 

The ratio TijTx depends on the angle between );„ 
and g cxt . For illustration, Fig. 5 shows its behaviour as a 
function of M e « for supersonic bodies moving along the ex- 
ternal gravitational field. In such a case, T2 is positive at 
large Mach numbers, implying that it contributes to drag 
the body. It becomes zero at M c s — 1.75 and negative for 
lower values. The relative contribution of T>2 becomes more 
important when approaching the transonic motion. At large 
Mach numbers, the ratio J-2/F1 increases monotonically but 
very slowly. 

In order to visualize the importance of including Fdf,2, 
Figure 6 shows the total drag force exerted on the perturber 
as a function of M, together with Fdf,i, in the deep-MOND 
limit (Lo = 1, To = 1/6). As anticipated, the sign of Fdf,2 



and g c 



In contrast to 



depends on the angle between v c 
the axisymmetric case, when and g axt are perpendicu- 
lar, Fdf.2 is positive for transonic Mach numbers and be- 
comes negative (reduces friction) at high Mach numbers. 
The contribution to the drag by the component T>2 is more 
important in the axisymmetric case but it is only notice- 
able (> 20%) at 1 < A4 < 1.5. Our numerical calculations 
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Figure 4. Distribution of the dimensionless integral I2 at cut-off planes y = (upper), x = (middle) and 2 = 10r p (lower) for 
= 0.53 (left column) and M = 1.25 (right column), in linear scale. The velocity of the perturber is along the 2-axis, whereas the 
external gravitational acceleration is along the x-axis and Lq = 1. 



show that Fdf,i scales with the size of the box domain as 
oc lnz max , whereas Fdf,i increases somewhat slower with 
Zmax- Therefore, the relative importance of Fof,2 is expected 
to be less for larger z max . 

Now, we wish to compare the drag force in deep MOND 
(1/0 = 1) and in Newtonian gravity (Lo = 0). Figure 6 also 
shows the Newtonian drag force experienced by a body of 
mass M, travelling at Mach number M when the wake has 
the same extent as in the MOND CEISG, 1.6. Zma,x — 
The drag force in MOND is a factor ctfj,g 2 larger than in 
Newton, where a is a form factor that depends on the Mach 



number and on the angle between v x and g cxt . In the ax- 
isymmetric case, a ~ 0.6 at low Mach numbers (M < 0.5), 
and becomes a ~ 0.5 at 0.5 < M < 1.0. For supersonic 
Mach numbers, a ~ 0.4 at 1 < M < 1.5 and increases 
monotonically with Mach number up to ~ 0.8 at high Mach 
numbers. The explanation for a — 0.8 at high Mach numbers 
is covered in detail in the Appendix C. When the direction 
of motion is perpendicular to the external field, a is very 
similar to its value in the axisymmetric case at M < 1. At 
supersonic velocities, a ~ 0.6 at 1 < A4 < 1.5, and falls 
monotonically down to 0.5 at high Mach numbers. Hence, 
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the drag force may vary with the angle between Voa and g oxt 
by as much as 50 percent. 



5 SOME IMPLICATIONS 

5.1 Wakes by galaxies and falling groups in 
clusters 

In this section, we discuss the implications of X-ray obser- 
vations of the morphology of wakes in clusters of galaxies 
for modified gravities. As in Furlanetto & Loeb (2002), let 
us assume that the galaxy or group of galaxies moves super- 
sonically through a constant density cluster core surrounded 
by a isothermal envelope: 



Pic(r) = 



Po 



if r < r c 



2p /[l + (r/r c ) 2 ] if r > r c 



where po is the density in the core and r c is the core radius. 
The emitted surface brightness S = J (e///47r)cH, where e// 
is the bremsstrahlung free-free volume emissivity, will be 
enhanced in the wake by: 



SS 



ST. 



(36) 



S ~ X E ' 

where E is the column density, and \ ~ 4 if the gas is 
isothermal, or w 5 if it varies adiabatically (Furlanetto & 
Loeb 2002). Maps of X-ray for a M = 1.25 body with v x 
perpendicular to the external gravitational field are shown 
in Fig. 7. The structure of the emission in the wake is very 
similar to that formed by a "compact" gravitational per- 
turber under Newton gravity. The projected X-ray emission 
along line-of-sights perpendicular to the direction of motion 
is roughly independent of the location in the wake, except 
near the edges of the cone, and is given by 



5S_ 
S 



2tt X 



GM 



(2 + tt) noPc OG T*c 



(37) 



As a consequence, in searching for the wake, one may expect 
an abrupt jump in surface brightness at the edge of its cone. 

We derived the gravitational wake induced by an ex- 
tended body with a Plummer profile having a fast density 
decay at large radii (p oc r -5 ) in an attempt to model the 
baryonic mass of a certain bound object. Wakes in the cold 
dark matter (CDM) scenario are expected to be different 
than in purely baryonic MOND because the dark matter 
component in the halo of galaxies decays as r~ 2 , much more 
slowly than the baryonic mass density (< r~ 3 ). To illustrate 
this, Figure 7 also shows the X-ray emission generated by a 
pseudo- isothermal perturber with core radius r v in the New- 
tonian case. We see that the X-ray emission in the standard 
CDM scenario is more cuspy. A MOND wake can be distin- 
guished from a CDM wake by detecting a sharp X-ray en- 
hancement along the Mach cone. Furlanetto & Loeb (2002) 
made a detailed analysis of the wake morphology in CDM 
models for collisionless and fully collisional (fluid) dark mat- 
ter (FDM) in the supersonic case. They found that, because 
in the collisional case the dark halo is truncated by ram pres- 
sure stripping, the X-ray emission of the wake is rather flat, 
similar to that we find in MOND. Due to the quantitative 
similarity between the wake in MOND and in FDM, many of 
the observational suggestions raised by Furlanetto & Loeb 
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Figure 5. Ratio between and T\, as defined in Eq. (35), 
versus the effective Mach number, for the axisymmetric case. At 
M e S < 1, the ratio depends on the adopted value for Lq and 
hence is not shown. 



to distinguish between FDM and CDM can be used to dis- 
tinguish between dark matter or MOND. We can repeat the 
reasoning of Furlanetto & Loeb (2002) and argue that the 
observations of the wake of the elliptical galaxy NGC 1404 
in the core of the Fornax group marginally support CDM 
against MOND, but the evidence is very weak. In the last 
decade, this type of observations has improved considerably 
(e.g., Drake et al. 2000; Neumann et al. 2001; Machacek et 
al. 2005, 2007; Sun et al. 2006). Still, it is difficult to draw 
some firm conclusions because of the difficulty to isolate the 
structure of the gravitational wake from the hydrodynamical 
wake, that is the mass in the wake stripped from the own 
galaxy by ram pressure. We must warn that MOND and 
FDM predict the same structure of the gravitational wake 
past a galaxy but, in many other astrophysical aspects, they 
must give different predictions since they are not equivalent. 

In principle, observations of the tails of subclusters of 
galaxies are a potential route to distinguish betweeen colli- 
sionless CDM and MOND. However, it is a well established 
issue that MOND still requires dark matter at cluster scales 
(e.g., The & White 1988; Sanders 1999). The inclusion of an 
isothermal dark matter component in MOND erases some- 
what the abovementioned differences between the MOND 
wake and the CDM wake. The observed displacement be- 
tween the X-ray peaks and the associated mass distribution, 
as derived from lensing data in the Bullet Cluster, basically 
rules out FDM (Markevitch et al. 2004) but not necessarily 
MOND (Angus et al. 2007). 



5.2 DF timescale in a spherical system 

In §4.2 we derived the DF force experienced by a body of 
mass M travelling on a rectilinear orbit through a homoge- 
neous fluid medium in deep-MOND and found that: 



Fdf 



a 4ttG 2 M 2 po 
Po v l° 



In A, 



(38) 
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Figure 6. DF force in a gaseous medium as a function of the 
Mach number when !)„ and g cxt are parallel (top) and orthog- 
onal (bottom). The open diamonds correspond to the drag by 
the component T>i in the deep-MOND limit and the asterisks 
correspond to the total drag. The dashed line and the triangles 
show the drag in Newtonian gravity. To help comparison, the to- 
tal drag in the axisymmctric case and the Newtonian drag have 
been also plotted in the lower panel (dotted line and dashed line, 
respectively) . 



where a depends on the Mach number and on the angle 
that makes the velocity of the perturber and the direction 
of the external field. How does it compare with the DF force 
in a collisionless medium? The formula for the collisionless 
MOND DF force was derived by Ciotti & Binney (2004). 
They show that it is similar to the Newtonian case but re- 
placing G — > Gao/3ext plus an extra factor of \f2: 



FDF,no-col 

where 
H = erf 



1 4ttG 2 M 2 po 



V2c 



2 Voo 

exp 



fflnA, 



2a 2 



(39) 



(40) 



with a the velocity dispersion of the Maxwellian distribu- 
tion of velocities of background particles (e.g., Binney & 
Tremaine 1987; Sanchez-Salcedo et al. 2006). By comparing 
Eqs (38) and (39), we can generalize the conclusion of Os- 
triker (1999) but for MOND: since the functional form of the 
gaseous DF drag is much more sharply peaked near A4 = 1, 
the drag force near A4 = 1 is larger in a gaseous medium 
than in a collisionless medium with the same density and 
a = Coo- At M w 1, using a « 0.5 (see §4.2) there is factor 
of ~ 3 difference in the MOND force between the two cases. 
For M < 1, the drag force is generally larger in a collision- 
less medium than in a gaseous medium, because pressure 
forces in a gaseous medium create a much more symmetric 
density perturbation in the background. 

For many problems of astrophysical interest, it is con- 
venient to have the DF inspiraling timescale for a body 
that is initially on a circular orbit. For circular orbits, the 
above formula is correct as long as the maximum impact 
parameter in the Coulomb logarithm is taken as ~ 2_R P , 
where R p is the instantaneous orbital radius of the perturber 
(Sanchez-Salcedo & Brandenburg 2001; Kim & Kim 2007). 
Consider a massive body embedded in the gaseous outer 
spherical envelope of a galaxy of mass Mg, with circular 
speed v c = (GMcao) 1 ' 4 . Suppose that the gas is isother- 
mal and in hydrostatic equilibrium in the gravitational po- 
tential of the parent galaxy. In a typical galaxy, the outer 
parts are expected to be in the deep-MOND regime, hence 
ffoxt = VGMGdo/r and no = (GMg/cio) 1 ^^ 1 . For the sake 
of clarity, let us assume that the gas has the virial temper- 
ature: Coo — i> c /\/2- Imposing hydrostatic equilibrium, the 
unperturbed gas in the envelope will pursue the following 
profile 



(41) 



where p a is the gas density at the radius of reference r a . 
By substituting the values for /io and p g into Eq. (38), 
and equating the rate of decrease of angular momentum 
d(Mv c r)/dt to the torque rF DF , we find the deep-MOND 
evolution of the orbital decay of a massive perturber's near- 
circular orbit in a gaseous isothermal spherical distribution 



r 

^"init 



exp 



t-tin 
Tm 



where n n i t is the orbital radius at t = t ln i t and 



tm = 



Ps 



1 



a 47rG 2 Mp a lnA' 



(42) 



(43) 
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with fi s = fJ,o(r s ). 

Consider now its equivalent Newtonian system (ENS) 
that is, the Newtonian system (dark matter plus gas) in 
which the baryonic gas has exactly the same density distri- 
bution as in the MOND system. The density distribution of 
dark matter in the ENS is: 



In the ENS, the dark matter component added to have the 
same "dynamics" also contributes to the DF experienced by 
the body. The frictional force in the ENS is 



Fdf 



4nG 2 M 2 (p 3 +0.428p d 



■In A. 



(45) 



The factor 0.428 arises because the dark matter has been 
considered collisionless (see, e.g., eq. [7-25] of Binney & 
Tremaine 1987). Equating again angular momentum loss 
with the torque, we find tens, the time that the body takes 
to reduce its distance a factor e in the ENS: 



0.43 



tens 



1 + 0.428^ 9 



r 2 ■ 

1 init 



4nG 2 M p a In A' 



(46) 



where lZ a 



Pdm/Pg 



=const. The ratio between the DF 



timescale in MOND and in ENS is: 

= 4.6(1 + 0A28TZ g )p 2 nit . (47) 

tens 

Here /ii n j t = ^to(nnit) and we used that a ~ 0.5 for M — 1.4 
(see §4.2). In terms of ^t^nit defined as the ratio between 
the mass in dark matter and the baryonic mass (stars plus 
gas) inside n n i t , we have 

1 



TM 
TENS 



= 4.6(1 + 0.428^ 9 



[1 + TZt 



(48) 
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In dwarf protogalaxies, before the gas is turned into stars, 
lZ g ~ IZt, and for typical values of dark matter contents 
in these systems (JZ t — 20), the gaseous DF timescale in 
MOND is 10 times shorter than in the ENS. Condensed 
objects that form early (e.g., globular clusters) could spiral 
into the centre of their host galaxy more rapidly than would 
be predicted in the standard CDM haloes. 



6 SUMMARY 

After a discussion about the Bondi-Hoyle problem, we have 
calculated the gravitational density wake of a perturber 
moving through a uniform gaseous medium in MOND. A 
hydrodynamical treatment provides useful insight into the 
problem of DF. In order to describe the response of the far- 
field medium to a small perturber, we focused on the case 
when the perturbation is dominated by a constant external 
field. The analytical knowledge of the wake structure in this 
simple case is useful to interpret fully non-linear simulations 
of the MOND dynamics of realistic systems. The structure 
of the wake depends on the angle between the velocity vector 
and the external gravitational field. We have considered two 
cases: the velocity of the perturber being either parallel or 
orthogonal to the external field. For an intermediate case, 
the structure of the wake will lay somewhere between the 
extremes described here. The density wake is decomposed 
into two linear contributions. The morphology of the domi- 
nant contribution is a scaled version of the Newtonian wake, 
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Figure 7. Normalized X-ray surface brightness maps for a 'com- 
pact' body moving at M = 1.25 in the deep-MOND non- 
axisymetric case for a line-of-sight along the j/-axis (upper panel) 
and along the ai-axis (middle panel), along with the X-ray map 
for an 'extended' pseudo-isothermal body in Newtonian gravity 
(bottom). The background field lays along the x-axis. The ob- 
server's line of sight makes an angle 7r/2 with respect the velocity 
of the perturber. 



being the MOND perturbed density a factor p,^ larger than 
in Newtonian gravity. The second contribution to the wake 
depends greatly on the angle which the direction of motion 
makes with the external field. 

The MOND DF force on the perturber induced by its 
own wake, as a function of the Mach number has been also 
derived and compared with the drag force in the Newtonian 
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case. It is important to know the dependence of the drag 
force strength on Mach number to study the circularization 
of orbits. Our results confirm earlier analyses suggesting that 
the DF force is higher in MOND than in Newtonian grav- 
ity. The DF drag is larger when the motion of the perturber 
is along the external field, especially at Mach numbers be- 
tween 1-1.5. In the context of the sinking satellite problem, 
the recent claim that the existence of an extended globular 
cluster population in Fornax is problematic for MOND gains 
strength (Sanchez-Salcedo et al. 2006; Nipoti et al. 2008). 
We find that MOND predicts an even faster orbital decay, 
especially when the satellite arrives to the halo of a galaxy 
before the gas has turned into stars. 

In MOND we show that a supersonic perturber gen- 
erates a wake in the gas with a well-defined Mach cone in 
which the the surface density increases substantially in a 
narrow region and then flattens. Because CDM haloes ex- 
tend to larger radii than the baryonic mass, the wakes in the 
two models have significant morphological differences. 

Inherent to the Chandrasekhar treatment of DF is the 
assumption that the perturber interacts with unbound parti- 
cles. In MOND, this requirement immediately demands the 
inclusion of the external field. Once the background gravi- 
tational field is included, two-body orbits of distant parti- 
cles are unbound, the interaction of the perturber with the 
medium is quasi-Newtonian, and the effects of distant en- 
counters can be simply added. This topic is closely related 
to studies of relaxation processes in self-gravitating media. 
In a forthcoming paper we will present a standard derivation 
of the DF drag in collisionless systems in MOND. 
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APPENDIX A: THE BONDI-HOYLE RADIUS 
IN MOND 

For a point-like particle, the linear MOND field equation 
(Eq. 4) is readily solved (Milgrom 1986): 

1 GM 



$i(r) 



Mo V(l + £<>)(x 2 + y 2 )+z 2 ' 



(Al) 



where we have supposed that the external field is along the 
2-axis. We are now in a position to derive the Bondi-Hoyle 
radius tbh,\ for the potential $i. We know, however, that 
<E>i is not a good approximation in the vicinity of the particle. 
Nevertheless, since |$i| overestimates the depth of the po- 
tential once we approach the perturber (i.e. |$i| > |$ OX act|, 
where <3>cxact is the exact potential for a point mass), the 
Bondi-Hoyle radius derived using the potential $i provides 
an upper limit to the exact Bondi-Hoyle radius. 

In the rest frame of the perturber, consider a stream- 
line that begins at x, y, z coordinates (x, y, — oo). A straight 
streamline is susceptible to bending and shocking if the ve- 
locity of the gas on that streamline (adding the thermal 
velocity) is larger than the maximum value of the escape 
velocity along that streamline: 

2 GM - 2 2 (A2) 



M 0(l + i ) 1/2 y^T 



>vL+c 00 



where we have assumed that Voo and g ext lie both along the 
z-axis. According to the above equation, all the streamlines 
that start with impact parameter b such as: 

2 GM 



b < 



(A3) 



bend significantly and go through shocks. In analogy with 
the Newtonian case, the MONDian Bondi-Hoyle radius is: 

1 GM 

If g cxt lies along the a>axis (so that g ext and Voc are 
perpendicular), the condition for blending is: 
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GM 



Mo ^x 2 + (1 + L )y 2 
The streamlines subject to terminate in a shock are 



(x 2 + (l + L )y 2 ) 1/2 < 



GM 



Ho (4(1 + X s 



(A5) 



(A6) 



We see that, due to the dilation along the axis of the back- 
ground field, the curvature of streamlines with the same 
impact parameter depends on their azimuthal angle in the 
(x, y) plane. Nevertheless, there is a guarantee that the re- 
sponse of the gas will be linear at any radius beyond Itbha 
where 



1 



GM 



Ho cL(l +M 2 



(A7) 



In the limit M >• 1 and in the steady state, see e.g., Ostriker 
(1999), 



Fdf,n oc In 



v 



oc \ ' mm 



(C3) 



Therefore, for z max S> r m i n , we have Fdf,n oc and 
2 



F D F,N{M c tt) 



) Fdf,n{M) 

M e B ) 

(1 + Lo)F df ,n(M). 



(C4) 



Substituting into Eq. (C2), we find that the drag force ex- 
erted by T>i in the supersonic axisymmetric MOND case, 
is about (1 — To)/L6q 2 = (5/6)/i,7 2 larger than in Newtonian 
gravity. 



APPENDIX B: THE COMPONENT V 1 IN THE 
AXISYMMETRIC CASE FOR A 
TIME-DEPENDENT PERTURBATION 

We may be interested in the temporal evolution of the wake 
when the perturber is placed at t = 0. The Newtonian case 
was studied by Ostriker (1999) who found that the induced 
wake density for this finite-time perturbation has the follow- 
ing analytic form 

w) = *£^k, (bi) 

sjz 2 -B 2 R 2 
where B 2 = M 2 - 1 and 



2 if m > i, ir + (z - vovty > (cooty 

z/R > 8, and z < fic^t/M; 
1 if R 2 + (z-v aD t) 2 < ( Coo t) 2 ; 
otherwise. 



We can take advantage of this Newtonian result to get T>\ (t) 
in MOND, at distances larger than a few r p , as: 



2>!(t) 



1 - T Q 

= ^fo^GM/^ 
sjz 2 - B 2 R? 

Here we used again that M\ = Mjr\. 



APPENDIX C: THE DRAG FORCE IN MOND 
AND NEWTON AT HIGH MACH NUMBERS 

According to Eq. (33), the drag force experienced by a Plum- 
mer body in Newtonian gravity is: 

Fdf,n = — 2 — T\ . (CI) 

^ c oo 

Let us compare Fdf,n with Fdf,i at M ^> 1 in the axisym- 
metric case. Note that at those Mach numbers Fdf,2 can be 
ignored. Combining Eqs (33) and (CI), the drag force com- 
ponent Fdf,i(M, Zmax) exerted on a body moving at M by 
a wake with extent 2 max can be recast: 

FDF,1 = — o7i — I — r — \FdF,N {M cB , Z m ax) • (C2) 



